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Recent Type la supernovae jj and the CMB data 0 
show that our universe is spatially flat (f= 0) and 
accelerating at present time. Specifically, the experimen¬ 
tal situation at present is such that the present energy 
density consists of approximately 1 /3 of ordinary matter 
and 2/3 of dark energy with a negative pressure. For this 
dark energy to accelerate the expansion, its equation of 
state w = p/p must satisfy w < —1/3. 

The most promising candidate for this cosmic acceler¬ 
ation and dark energy seems to be a “true” cosmological 
constant (CC) A (with A > 0), for which the equation of 
state is simply w = — 1. Another interpretation of dark 
energy relies on the existence of an ultralight scalar field 
4> (“quintessence”) 0, in which case the dark energy den¬ 
sity is the result of the scalar field slowly evolving along 
its effective potential V (</>), with — 1. Some alternate 
explanations for the present acceleration beyond dark en¬ 
ergy include theories in more than four dimensions pj as 
well as the Chaplygin gas fE]. 

If the presently dominating dark energy and the ac¬ 
celeration of the universe are due to the “true” CC or 
“quintessence”, the expansion will keep accelerating per¬ 
petually, rapidly approaching a de Sitter (dS) regime. 
This implies a fate for open (f2*, > 0) and flat (f Ik = 0) 
universes in which the scale factor expands exponentially, 
a ~ e Ht , for an indefinitely long time. The same is des¬ 
tined even for a closed universe, unless 12 a is too low so 
that the universe collapses before the CC becomes pre¬ 
dominant. More importantly, an eternal expansion and 
acceleration would exhibit future event horizons, which 
presents a rather serious challenge to critical string the¬ 
ory S' Namely, the finite event horizon implies the ex¬ 
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istence of regions in space that are inaccessible to light 
probes, thereby preventing the definition of asymptotic 
states and therefore the construction of the conventional 
^-matrix description required for field theory. This im¬ 
plies the impossibility of formulating the S'-matrix for 
strings, which are by definition theories of the on-shell S'- 
matrix, in such backgrounds. A noncritical string theory, 
characterized by a “graceful exit” from the inflationary 
phase, has recently been considered in 0. 

Nonetheless, there are other models describing dark en¬ 
ergy in which the dS phase is transient, thereby bringing 
back compatibility with string theory. One loophole pro¬ 
viding an alternative to a perpetually accelerating uni¬ 
verse is to have a potential V(<j>) which has a minimum 
at </>o, with V(<t> o) = 0 and V (</>o) > 0. In this case, 
during the coherent oscillation phase, the energy density 
decreases as a~ 3 (w = 0), leading to a decelerating uni¬ 
verse. In another scenario, V((j>) has a minimum V((f> o) 
= 0 and then becomes flat for <fi > <fi o- After passing 
the minimum, the epoch of kinetic energy domination 
sets in where the energy density scales as a -6 (w = 1), 
resulting again in a decelerating universe. Nonetheless, 
the simplest case is to have a negative CC. In this case, 
no matter the universe is open, flat, or closed, or how 
extremely tiny the CC is, it collapses eventually, thereby 
bypassing the troublesome problems with the event hori¬ 
zon. As a way of example, we mention a mechanism to 
generate a negative CC in 0 (for supergravity theories, 
see, e.g., 0 ). 

In the present paper, we expand our previous work 
M on the RG running of the CC in the standard model, 
to a generic field theory model. We then study the fu¬ 
ture of the universe in such a generic scenario. Recent 
works narnia show that even a “true” CC cannot 
be fixed to any definite constant (including zero), ow¬ 
ing to RG running effects. Particle contributions to the 
RG running of A due to vacuum fluctuations of massive 
fields have been properly derived in M , with a some- 
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what unexpected outcome that more massive fields do 
play a dominant role in the running at any scale. If our 
chosen value for the RG scale p is always below the low¬ 
est mass in the underlying quantum field theory, then 
the scaling evolution of the CC A (p) towards p —> 0 is 
essentially given only by two parameters which depend 
on the field content (particle masses as well as the to¬ 
tal number of massive degrees of freedom including the 
+/- sign for bosons/fermions). The same structure of 
the scaling behavior is also expected for A arising from 
vacuum condensates of scalar fields 0 , in which case 
there is an additional dependence of the parameters on 
the particle interaction. We show that even if we start 
with a positive A in a flat universe, the scaling behavior 
may naturally be such as to provide a negative A for an 
indefinitely long time, thereby changing the ultimate fate 
of the universe. 

We have already mentioned the possibility that a flat 
universe may collapse in the future. In most models of 
dark energy based on extended supergravity, this occurs 
in the nearby future, within the next few billion years 
0 Although provided with a negative CC for an indef¬ 
initely long time, our scenario allows quite a distinctive 
feature: the universe never collapses, starting to deceler¬ 
ate relatively soon and increasing the scale factor to its 
maximum value in infinite time. 

In the formalism of quantum field theory in curved 
space-time, the cosmological constant receives divergent 
contributions originating from zero-point energies of the 
quantum fields considered. Being one of the parameters 
in the Lagrangian, the cosmological constant undergoes 
renormalization (for earlier works on the subject see e.g. 
01. After the appropriate choice of the renormaliza¬ 
tion scheme which correctly comprises the effects of mass 
thresholds 0 , one can formulate the renormalization- 
group equation in the form 
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where p is the renormalization scale, p\ represents the 
energy density attributed to the cosmological constant, 
and i is an index of any massive degree of freedom (cor¬ 
responding to the quantum field), both fermionic and 
bosonic. The quantity Ii acquires the value +1 for 
bosonic and —1 for fermionic degrees of freedom. It is 
clear that more massive fields dominate 0 at any scale 
and this departure from the Appelquist-Carazzone the¬ 
orem 0 results from the (mass) 1 * * 4 dimensionality of p\ 
iljj. The following comment is in order. It is the high 
dimensionality of the CC energy density that causes the 
failure of the decoupling theorem; the peculiar mass- 
dependent renormalization scheme we used in 0 in or¬ 
der to renormalize the zero-point energy agrees with the 
MS scheme in the high-energy limit. The result is fur¬ 
thermore supported by the failure of the decoupling the¬ 
orem in the case of vacuum condensate of the Higgs field 


0 where the mass-dependent renormalization scheme 
is unambiguous . 

In the approach described above, the particle zero- 
point energies are taken as a generic contribution to the 
running of the cosmological constant. The results for the 
running are obtained at the one-loop level. In the frame¬ 
work we used, quantum field theory of matter (particles) 
is defined on the curved classical background which is 
described through the vacuum action 

S = j - 2A ) + 

J d 4 X\/—g + a^R 2 + a 4 df?]( 2 ) 

This equation, however, does not include the possible 
contributions of quantum gravity. As discussed in djJ, 
one could safely neglect the gravity effects related to all 
terms in 0 proportional to a*, and the nonminimal term 
at the present time and scale. However, nonper- 
turbative effects of quantum gravity could play a certain 
role even today. 

Althogh it is presently well known that there are 
tremendous difficulties in the formulation of quantum 
field theory of gravity, it is still possible to renormal¬ 
ize the nonrenormalizable theories in the framework of 
effective field theories 0 A successful example of such 
a theory is a chiral perturbation theory as a low-energy 
nonlinear realization of QCD. Gravity is potentially even 
a better effective theory, since quantum corrections ap¬ 
pear to be rather small, persisting so up to the Planck 
scale. It is then possible, using the effective field theory 
to calculate the effects of quantum gravity which emerge 
from the low-energy part of the theory without knowing 
the ’’true” quantum field theory of gravity 0. 

What really makes the effective field theory of grav¬ 
ity different from other effective theories is the eventual 
existence of singularity 0 in the future - which, in the 
limit of extremely low energy (corresponding to the the 
wavelength probed of the order of the size of the uni¬ 
verse), breaks the validity of the theory in the infrared 
limit. For example, conformal invariance breaking (which 
is also broken by grav itons!) causes infrared instability 
in de Sitter space |2(| ■ Therefore quantum gravity may 
lead to strong renormalization effects at large distances 
owing to the infrared divergences 21]. 

Recently, Bonanno and Reuter [22] derived RG equa¬ 
tions with an infrared cutoff which at a given scale stops 
the running of A in the infrared. The basic framework is 
the effective average action T [g^u], which is a Wilsonian 
coarse grained free energy. It depends on a momentum 


1 Recently, Gorbar and Shapiro O studied the form of decou¬ 

pling for the CC energy density in quantum field theory in the 

external gravitational field on a flat background. The results fail 

to reproduce a correct high-energy behavior as given in the MS 
scheme, and, therefore, cannot be conclusive. 
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scale k, and defines an effective field theory appropriate 
for the scale k. An additional hypothesis in the formula¬ 
tion of the RG behavior of A and G is that, in the limit 
k —> 0, both A (k)/k 2 and k 2 G{k ) run into an infrared 
attractive non-Gaussian fixed point. 

The scale k plays a role of an infrared cutoff 0]. Intu¬ 
itively, an infrared cutoff is necessary since, for example, 
at the present time t = t 0 (t = 0 at the Big Bang), one 
has to integrate quantum fluctuations down to the scale 
k > A; larger scales for which k < A should obviously 
be truncated. 

However, it is very difficult to imagine the precise phys¬ 
ical mechanism which acts as a cutoff. It would depend 
on the problem under consideration and may include 
many scales including, e. g., particle momenta, field 
strengths, the curvature of spacetime, etc. We address 
this problem and our choice of the scale in the discussion 
later on. 

The integration of Eq. 0 gives a closed form of the 
cosmological constant energy density at any value of the 
renormalization scale p: 

PaM-Pa(O) = E JiTO ^ln(l + ^2) • ( 3 ) 

At renormalization scales much smaller than the mass of 
the lightest quantum field the expansion in powers of /i 2 
gives the following form for p\\ 


Pa(m) - Pa(0) 
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In the equation given above, Nf, and Nf denote the to¬ 
tal numbers of bosonic and fermionic massive degrees of 
freedom, respectively. 

The problem to be solved next is the choice of the 
renormalization scale p in cosmological settings. In 
an evolving universe, it would be natural to identify 
the renormalization scale p, representing typical parti¬ 
cle (graviton) momenta in the universe, with some en¬ 
ergy scale characteristic of the stage of evolution of the 
universe. During epochs of radiation domination, associ¬ 
ation of the scale p with temperature, p ~ T, is by all 
means straightforward. In our approach, we adopt as a 
natural choice 


p = p 1 '\ (5) 

where p stands for the total energy density of the uni¬ 
verse. This choice determines the renormalization scale 
in terms of a quantity well defined in cosmological set¬ 
tings. Using 0 we obtain a consistent definition of the 
scale /i for all epochs of the evolution of the universe. 


Moreover, in the radiation dominated universe, the defi¬ 
nition of the renormalization scale 0 coincides with the 
straightforward association p ~ T mentioned above. 

Next we discuss our choice of the scale, p = p 1 / 4 , and 
compare it with other possible choices, especially with 
that proposed by Shapiro and Sola [uj, p ~ VR ~ H , 
where R is the curvature scalar and H denotes the Hubble 
parameter. 

Let us recall that in the decoupling theorem, it is the 
momentum of particles entering a quantum loop, and the 
mass of particles inside the loop, that are to be compared 
(with massless gravitons on external legs in the case of 
CC). On the other hand, the coupling constants play no 
role in the decoupling. It has been argued 0 that the 
value of the curvature scalar R can be taken as an order- 
parameter for the gravitational energy, i.e., p ~ y/~R. The 
Einstein equations imply that, at present, 

Po ~ Rl /2 ~ (Gpo) 1/2 ~ Ho , (6) 

where G denotes the Newton’s constant. According to 
the choice 0 . the relevant scale cannot be less than 
that of typical momenta of background gravitons, which, 
at present, are of the order of 10~ 4 eV and not of 
H 0 ~ 10“ 33 eV. We notice that the large separation 
of the scales 0 and 0 is actually due to a coupling 
of gravitons to ordinary matter (~ GE 2 with E being 
a typical graviton energy), inherently present in the ex¬ 
pression 0 . Indeed, one finds that the ratio of scales 
© and © is essentially given by the square root of the 
graviton coupling constant. It is thus pretty obvious that 
the Shapiro-Sola choice (5) means a restriction to “soft” 
graviton momenta, whereas our choice (4) involves typi¬ 
cal (or “hard”) graviton momenta. 2 In turn, this implies 
that the choice © corresponds to CC’s that are nearly 
uniform in both space and time, thus being compatible 
with models which treat A as an (almost) true constant. 
In this case, A acts as a source of gravitons with ex¬ 
tremely large wavelengths, of the order Hq . Moreover, 
with the choice ©, the running of the CC on cosmolog¬ 
ical scales is essentially negligible. On the other hand, 
adopting the choice ©, one arrives at a model with the 
pronounced running of A, comparable with that of ordi¬ 
nary matter. This feature is clearly visible in our figures, 
see below. 

Let us in more detail discuss our advocacy for the scale 
as given by 0 . First of all, we would like to have the 
correct limit towards the past radiation epoch, where all 
relevant particle scales are essentially given in terms of 
the temperature of interacting particle species. Obvi¬ 
ously, this is not possible with the choice ©. In addition, 
in our approach the energy-momentum tensor is no 


2 The above separation of the scales resembles a distinction be¬ 

tween hard (~ T) and soft (~ gT) momenta encountered in 

quantum field theories at finite temperature/density, where g is 
the gauge coupling constant of the underlying theory l23l . 
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longer separately covariantly conserved, but only is the 
sum T M j, + (see (UTl) below). This implies that 

there exists an interaction between matter and the CC 
which causes a continuous transfer of energy from matter 
to the CC or vice versa, depending on the sign of the in¬ 
teraction term. We shall see below that in our approach 
the CC and matter are always tightly coupled, with a 
consequence that the CC actually behaves more like ordi¬ 
nary matter than a “true” CC. Consequently, both com¬ 
ponents are assumed to act as short wavelength sources 
of gravitons with p » H ~ 1/t ( H ~ 1/fis adequate 
only for the CC with a negligible interaction with ordi¬ 
nary matter). A similar distinction of scales can be found 
in the paper [2^, with the statement that the choice 
p ~ H ~ 1/t is adequate only for the case of perfect ho¬ 
mogeneity and isotropy, otherwise new scales much larger 
than the cosmological one need be introduced. Finally, 
the choice p ~ H ~ 1/t, in our approach, which aims 
to take only effects of matter on the renormalization of 
A, causes a few problems on their own. For instance, if 
one identifies the parameter A (CD and 0 below) with 
the square of the highest mass in the underlying particle 
theory, one can fit the present experimental data for pa, 
Pa(p = po) ~ 10 -12 eV 4 , only with A ~ AfJ> ; , which is by 
no means plausible(obviously, the Planck scale is natural 
in the context of quantum gravity but not in the context 
of any particle theory). The only remedy would be to 
give up the “nice” infrared limit, pa(0) = 0, which ne¬ 
cessitates antropic considerations to explain the nonzero 
infrared value. Furthermore, as already stated, in this 
case one finds that p(p = po) ~ pa(0), which leads to 
totally negligible running, therefore rendering this case 
completely uninteresting for the purpose of the paper, 
i.e., for the study of the fate of the universe. 

Clearly, the choice 0 constitutes a viable model as 
long as the total energy density is positive. Therefore, we 
consider the evolution of the universe in a time interval 
in which the identification of Eq. 0 is sensible. As we 
shall see below, the corresponding time interval turns out 
to be the whole interval to < t < oo, and therefore the 
choice of the renormalization scale given by 0 actually 
poses no restriction at all. 

From this point on we consider the evolution of the 
universe in cosmology with the running (indirectly time- 
dependent) cosmological constant energy density. We 
assume that the present value of the cosmological scale 
po ~ 2 x 10~ 3 eV is sufficiently smaller than the lightest 
of the masses of the fields in the theory. In this case, 
the asymptotic expansion © applies. In view of recent 
observations of neutrino oscillations by Superkamiokande 
111 and SNO [25], we find the expansion 0 justifiable 
as these experiments seem to indicate that neutrinos do 
have a significant mass above 10~ 3 eV. In addition, we 
assume that the infrared limit of the cosmological con¬ 
stant energy density vanishes, pa(0) = 0, hence avoid¬ 
ing anthropic considerations to explain this value. In 
this way, we arrive at a model in which the cosmolog¬ 
ical constant energy density is dynamically determined 


at any scale. Although the coefficient of the p 2 -term 
in (3) is severely restricted by the big bang nucleosyn¬ 
thesis results [id . this term, if nonzero, still dominates 
at sufficiently low scales. Therefore, at sufficiently low 
scales (po > p > 0), relevant to the study of the ul¬ 
timate fate of the universe (i.e., for the corresponding 
times to < t < oo), we describe the running of the cos¬ 
mological constant energy density by the first two terms 
in the expansion 0 

p A = Ap 2 + Bp A , (7) 

which together with the definition of the cosmological 
energy scale 0 gives 

p A = Ap 1/2 + Bp. (8) 

We have shown in the previous paper 10] that also for 
the vacuum-condensate contribution to A one should ex¬ 
pect the same scaling behavior of pa as in 0. Hence, one 
may consider 0 a scaling behavior for generic A, where, 
besides the definition in 0, the coefficients A and B now 
take on an extra dependence on the particle interaction. 
In order to round up the survey of our cosmological model 
there remains to identify other components of the total 
energy density apart from the cosmological constant en¬ 
ergy density. In a general manner, we express the total 
energy density and pressure of the content of the universe 
in a form 

P = Pa + Pm j 

P = PA+Pm, (9) 

where we have used the symbols p m and p m to specify 
the density and pressure of all other components of the 
universe, apart from the dark energy ones. Thus we as¬ 
sume that the running cosmological constant is the only 
source of the present acceleration of the universe, i.e., the 
only dark energy component. 

In order to couple the cosmological term (being a dy¬ 
namical quantity) with the Einstein field equations 

Gfu, + A= —87t GTp V , (10) 

we completely adopt the prescription outlined in a re¬ 
view paper [2g| to move A to the right-hand side of Uni- 
interpreting it, at the same time, as a part of the mat¬ 
ter content of the universe, rather than a pure geomet¬ 
rical entity. After doing so 3 , it is the effective energy- 
momentum tensor that has the perfect fluid 


! If one rather prefers to keep a variable A-term on the left-hand 
side of rroi . then introduction of some additional terms into im 
is necessary, which are functions of the scalar in a tensor-scalar 
theory. 
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form and therefore satisfies energy conservation: 

a ~r~ + 3(p m + p m ) = 0 . (11) 

da 

The equation of state obeyed by the cosmological term 
is then of the familiar form, p A = —p\, and that obeyed 
by ordinary matter may have a general form p m = wp m 
( although in our analysis we restrict ourselves to the 
pressure-free Friedmann models with w = 0). Note from 
CD that p m scales in a familiar way, p m ~ a , when 
A = constant as well as w = 0. 

Furthermore, we introduce a condition (at a(0) = ao) 
coming from present observational results on the share 
of dark energy in the total energy density at the present 
epoch: 

PA (a = a 0 ) = ^aPo , (12) 

with ~ 0.7. The condition given above determines 
the parameter A in © 

A = (n° A - B)p 1 0 / 2 . (13) 

With substitutions u p = p/ po and y = ln(a/ao), Eq. itTTl) 
acquires the form 



FIG. 1: The dependence of the energy density components 
of the universe (in units of po) on the scale parameter a for 
the parameter values B = 0.6, fl® = 0.7, and w = 0. For 
large values of the scale factor, p m tends to zero, whereas pa 
approaches p and tends to a constant nonzero value. The 
hierarchy p m < pa < p is maintained for all values of the 
scale factor. 


^ + 3(1 + w)((l - B)u p - (n° A - B)ul /2 ) = 0 , (14) 

while the initial condition is u p {y = 0) = 1. This equa¬ 
tion has the following solution: 


Up = — = 

Po 


1—§(1 — ^ aX 1 + w ) ln (— 
z \ao 


1 2 


■ ( 18 ) 


The expressions for energy densities pa and p m in the 
case B = 1 are 




Po 


n° A - B 1 -n° A ( a 


1 — B 1 — B \ao 




The expressions for energy densities pa and p m are 


(15) 


PA — PO 


0° — x(l + ^a)(1 ~ ^a)(1 + w ) l n ( — 
2 V a 0 


9 - ~0\2n , „,a 2 ( ln I ff_ 


+ ^l-nmi + w) 


a 0 


(19) 


PA — P0 


rtk-B 


1 - B 

(i + B)(n° A -B){i-n 0 A ) 


(i -sy 
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-§(l-B)(l+«,) 
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1 -n° A \ 2 1 ' a \- 3{1 - B){1+w) ' 


1 -B 
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(16) 


Pm — P0 


(nl-B)( i-nfr) /_« 

(i -b) V«o 


-§(1-B)(1+ W ) 


1 - B 


a o 


-3(i-B)(i+ w y 


(17) 


Pm — Po(l « A ) 


1- f(l - ^a)( 1 + w)ln (— 

2 \ao 


( 20 ) 


The definition of the cosmological energy scale requires 
that uj 2 > 0. This requirement for values B > fl A leads 
to the existence of the maximal value of the scale factor. 
In general, the expression for the largest value of the scale 
factor a m ax is 


Q>max — ^0 


B-n ° \ 30-BX1 + U,) 


( 21 ) 


while for the special case B = 1, we find 


In the special case B = 1, the solution of m has the 
form 


a — „ nP 3(1 -« a )(1 +“’> 
Umax — LLQ& ^ 


( 22 ) 





























6 



ln(a/a 0 ) 


FIG. 2: The dependence of the energy density components 
of the universe (in units of po) for the parameter values B = 
Ha =0.7 and w = 0. For large values of the scale factor a , all 
energy densities approach zero. The hierarchy p m < pa < p 
is valid for all values of the scale factor. 



FIG. 3: The dependence of the energy density components (in 
units of po) of the universe for the parameter values B = 0.85, 
Ha =0.7, and w = 0. The scale factor a has a maximal value 
a,max . For scale-factor values close to a mo .x, the pa becomes 
negative, while p m surpasses p. 



FIG. 4: The dependence of the energy density components of 
the universe (in units of po) for the parameter values B = 1, 
Ha = 0.7, and w = 0. The scale factor a has a maximal value 
a,max■ For scale-factor values close to a m ax, the pA becomes 
negative, while p m surpasses p. 



FIG. 5: The dependence of the energy density components of 
the universe (in units of po) for the parameter values B = 1.2, 
Ha = 0.7, and w = 0. The scale factor a has a maximal value 
dmax■ For scale-factor values close to a m ax, the pA becomes 
negative, while p m surpasses p. 


Finally, we can incorporate the laws of change of the 
total energy density with the scale factor into the Fried¬ 
mann equation for the flat universe (k = 0): 


8t tG 


-P- 


(23) 


Solving the Friedmann equation, we obtain the law of the 
evolution of the scale factor 


a = a o 




K-b 


K-B 


3(l-B)(l + m) 


For B = Ha, the a(t) law has the form 


( 24 ) 
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FIG. 6: Evolution of the scale factor a with time for five 
different values of the parameter B and = 0.7, w = 0. 
Notice that for values B < the universe expands without 
limits, whereas for B > the scale factor asymptotically 
tends to a m ax and reaches it at infinity. 



FIG. 7: The event horizon for the parameter values B < 0° 
in units 1/Ho. For values B > f1° , the event horizon becomes 
infinite. The parameter values used are 17° = 0.7 and w = 0. 
The vertical line at B = fl® is introduced as a guidance to 
the eye. 


1 + — Ho(t — f 0 )(l + «0(1 — fii) 


a = a o 

while for B = 1, we obtain 


3(l-n» )(l + m) 


(25) 


a = ao e 3(1 - S2 AK 1 +*") 




(26) 


An important feature of the evolution for B > is 
that, although there exists a maximal value of the scale 
factor, it takes the universe infinitely long to reach it. 
Moreover, for any finite time the total energy density p 
is nonnegative which makes our definition of the cosmo¬ 
logical energy scale generally acceptable (and not only 
for some finite interval of time starting now). 

The cosmological model described above also has in¬ 
teresting consequences for the event horizon dE{to) = 
n r°° dt' 

0 Jt 0 a(t') ■ 

In the case B = we can give the expression for the 
event horizon in the closed form 


2 1 

dE ^ = 2 — 3(1 — 52° )(1 + w)) ~Hq ' ^ 

Numerical integration for the values of the parameter 
B < gives the finite values for the event horizon as 
depicted in Fig. Q On the other hand, the event horizon 
for the parameter values B > diverges. This can be 
easily seen from the fact that in this parameter range the 
scale factor is bounded from below and from above: 


a 0 < a(t) < a max , (28) 

which leads to 

r°° dt' , , , r dt' 

act / — > ds(t o) >a 0 -. (29) 

J to a 0 J to a max 

As both the upper and the lower bound in the above 
inequality diverge, the event horizon is infinite in this 
parameter range. 

Let us now discuss the main features of our scenario 
in a more detail. As shown in our figures m several 
different outcomes concerning the fate of the universe as 
well as the scale factor evolution are possible, depending, 
in essence, on how the parameter B is related to f2°. 
Once B is picked, the parameter A is given by d. Our 
treatment of A and B as free parameters means that 
our analysis is, on the whole, model independent (i.e., 
we do not restrict ourselves to a given particle theory). 
Let us mention, for instance, that considering zero-point 
energies only, we have B = x i n standard model 
(treating neutrinos as massless). 

First, consider case B < as depicted in fig. |T| We 

see that p\ approaches a positive constant po ( 

when a —> oo, whereas p m approaches zero in the same 
limit. This means that the universe once entering the dS 
regime, stays there for an indefinitely long time, see fig. 
0 In this case, the running of p\ plays no essential role, 
since it is the asymptotic behavior of p\ and p m that 
matters in determining the ultimate fate. Of course, the 
























event horizon is finite, see fig. 0 Next, consider the lim¬ 
iting case B = f2° depicted in fig HI Now both pa and 
p m tend to zero when a —> oo. The hierarchy p m < Pa 
preserved in the whole running is crucial for the ultimate 
fate, which turns out to be the same as in the preceeding 
case, see fig. 0 The event horizon is also finite, see fig HI 
Now we consider the most interesting case B > as de¬ 
picted in hgsEE In all these cases, the universe left the 
dS phase relatively soon, when a ~ ooe 0 5-15 (depending 
on B ), then starting to decelerate. The reason for this is 
the change in hierarchy from pa > p m to p m > p\. Soon 
afterwards, when a ~ aoe 1-2 ' 5 (depending on B ), pa be¬ 
comes negative, staying such forever. We recall that for 
the “true” A (A < 0), the universe recollapses unavoid¬ 
ably. The reason why we have a different fate here is 
that the established hierarchy, p m >| p\ |, is preserved 
also asymptotically, where both components tend to zero 
energy density. Stated differently, since | pa | goes faster 
to zero than p m when a —> a max , the former component 
never has a chance to start dominating, a feature crucial 
for having recollapse within a finite time interval. Finally, 
we need to explain why a max < oo in this case, see fig. 
0 Stating formally, the recollapse here occurs at infinity 
(i.e., for t = oo), when | p\ | finally reaches p m . Since 
any recollapse (starting whenever) requires a ma x < oo, 
we have a(t = oo) = a max < oo, as seen from 1211221) 
and fig. 0 The matter dominance in the asymptotic 
regime, helped with the influence of negative A, makes 
the event horizon diverge, thus alleviating the problems 
with strings, see fig. 0 

The results on the evolution of the universe presented 
in this paper refer to the case of a flat universe (k=0). For 
this case, it is possible to obtain analytical expressions for 
the time evolution of the scale factor a(t). In the cases of 
open (k=-l) or closed (k=+l) universes, the curvature 
effects may become important in the considerations of 
the asymptotic evolution of the scale factor. In these 
cases, however, it is not possible to maintain the same 
level of analytical tractability. 

Finally, one would like to know if the above scenario 
is capable of solving the well-known problems that one 
encounters when explaining dark energy with a “true” 
CC. The first problem is to ensure a phase of inflation, 
an epoch when vacuum energy dominated other forms 
of energy density, characterized by p vac p°- Let us 
also mention the so-called coincidence problem, which 
requires an extremely large hierarchy between pf j and 
Pa at the epoch just after the inflation (contrary to the 
expectation from equipartition), in order to correctly re¬ 
produce pa today, with A = constant. Since both of the 
problems are addressed by “quintessence”, one may wish 
to know whether the above scenario could mimic some of 
quintessence models. This amounts to study the running 
© backward into a distant past. Whereas at extremely 
low energies we could keep only the first two terms in 
the expansion of © and perform the analysis without 
referring to any particular model, this is no longer pos¬ 
sible at high energies. Now we need all peculiarities of 


the underlying particle theory: particle masses, all rele¬ 
vant interactions, symmetry breaking properties, etc. For 
instance, the description of the running of the cosmo¬ 
logical constant energy density at high scales p requires 
the knowledge of the number and individual masses of 
the heaviest degrees of freedom (and not just a small 
number of specific combinations of these parameters like 
coefficients A and B). Such complete information is cur¬ 
rently not available even for the simplest theories, like 
the standard model. There is an additional technical 
complication in the region of high scales. Namely, one 
has to take into account the contributions coming from 
mass scales differing many orders of magnitude, which 
hinders a straightforward numerical treatment. In addi¬ 
tion, higher-loop contributions together with the running 
of masses and couplings, all of which we have obviously 
disregarded in ©, may no longer be irrelevant at high 
energies. We feel, though, that it would certainly be of 
importance to address the coincidence problem even in a 
toy model, in view of the obvious connection between p\ 
and p m as expressed in ©. Since our examples show a 
clear tendency for p\ to grow faster than p m for t < to, 
one may object that there is no room in our scenario to 
explain the discovery of the supernova SN 1997ff, which 
is believed to signal a switch from cosmological accelera¬ 
tion to deceleration [27j ■ Once again we must stress that 
this potential drawback could result from our ignorance 
of the full theory. We recall m that the smallness of the 
parameter A requires a fine-tuned relation among parti¬ 
cle masses of the underlying theory. One may therefore 
imagine that the running of masses, although believed 
to be weak at such energies and therefore ignored in (1), 
could easily change the sign of A , thereby depleting p\ 
at higher energies. 

As a last step, it would be very useful to compare the 
effects originating from the matter sector (treated in this 
paper) with quantum gravity effects on the running of 
the cosmological constant. Firstly, we could interpret 
the assumption pa( 0) = 0 as an effective inclusion of 
the effects of quantum gravity in our model since within 
our model there is no mechanism or principle that would 
dictate vanishing of p\ at p = 0. This choice certainly 
represents a rather crude approximation of the effects of 
quantum gravity, but, in our view, at the same time sin¬ 
gles out the most important of these effects. Secondly, in 
the absence of complete understanding of quantum grav¬ 
ity effects, especially in the infrared region of the scale 
p, it is not possible to say which type of effect (matter 
sector or quantum gravity) would dominate the running 
of the cosmological constant energy density at infrared 
scales. Still, our findings of significant (and possibly sign¬ 
changing) running of the cosmological constant energy 
density coming from the matter sector make reasonable 
the assumption that quantum gravity effects (apart from 
those setting pa(0) = 0) do not significantly change the 
running of the CC and the respective consequences on 
the evolution of the universe. We see our model as com¬ 
plementary to quantum gravity approaches, especially to 




9 


those which draw conclusions from the pure gravity the¬ 
ories assuming that the matter sector does not alter the 
pattern from the pure gravity theory [2^ |. Our results 
show that one cannot a priori expect the effects of the 
matter sector to play only a secondary role. We expect 
that only a full solution of quantum effects in gravity 
along with a careful treatment of contributions coming 
from the matter sector will resolve the problem of the 
running of the cosmological constant unambiguously. 

In conclusion, we have studied the renormalization- 
group running of the cosmological constant at scales well 
below the lowest mass in a generic particle theory. Hav¬ 
ing specified the renormalization-group scale in cosmo¬ 
logical settings, we considered the ultimate fate of the 
universe in such a scenario. As its main feature, our sce¬ 
nario naturally provides the possibility that the universe 
endowed with a positive cosmological constant nowadays, 


necessarily exits a de Sitter regime a few Hubble times 
afterwards. More surprisingly, although such a universe 
exhibits the maximium value of the scale factor, it never 
recollapses. Although the running into past times cannot 
be provided in a model-independent manner, the relation 
between p\ and p m implied in m indicates that, treated 
in a specific particle physics model, our approach might 
also contribute to the resolution of the coincidence prob¬ 
lem. 
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